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FOREWORD 


This report is No. 8 of a series written for the Diffusion 
Project presently being conducted at Colorado Agricultural and 
Mechanical College for the Office of Naval Research under 
Contract N 9 onr 82401. The experimental phase of this project 
is being carried out in a wind-tunnel at the Fluid Mechanics 
Laboratory of the College, The project is under the general 
supervision of Dr. M. L,. Albertson, Head of Fluid Mechanics 
Research of the Civil Engineering Department. 

To Dr. M. L. Albertson, and to Dr. D. F. Peterson, Head of 
the Civil Engineering Department and Chief of the Civil Engineer- 
ing Section of the Experiment Station, as well as to Professor 
T. H. Evans, Dean of the Engineering School and Chairman of the 
Engineering Division of the Experiment Station, the writer wants 
to express his appreciation for their kind interest in the pre- 
sent work, 

The writer also wishes to thank the Multigraph Office of 


the College for the able service it has rendered. 
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On the Asymptotic Behavior of any Fundamental Solution of 
the Equation of Atmospheric Diffusion and on a 
Particular Diffusion Problem 
by 
Chia-Shun Yih 


Abstract 
In this eanedy the asympotic behavior of any fundamental 
solution of the differentaal equation of atmospheric diffusion 
is studied, It is found that if the wind velocity and the 
diffusivity increases monotonically with height, then the 
"amplitude" and the spacing of the zeros of the fundamental 
solution will decrease asymptotically in certain définite 
ways. AS an application a particular problem in atmospheric 
diffusion is solved at the end. 
1. Introduction 
If one neglects the longitudinal diffusivity in comparison 
with the vertical ieee uel GLE the ecuation of diffusion can be 


written as ean 
=" Gl! Pas slitz e C. \ 
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where c, U, and K are respectively the concentration, the 
wind velocity and the vertical diffusivity. =; and K being 
functions of y, only, and x) and Yi being measured 
respectively in the horizontal and the vertical directions, 
With ¢, denoting the ambient sonnentnaticn, and h 
denoting a reference length, the last equation can be written 


in the dimensionless form: 
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C=C, 
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Q 
U. being a certain reference velocity. In sean hea ob 


2) 
fusion, u and D are usually assumed to be monotonically 
increasing functions of y. 

To solve the differential system consisting of (1) and 
various boundary conditions, the method of separation of vari- 
ables will be used. Assuming 


?= Xa Yon eu 


and substituting in (1), one has 
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where the primes denote differentiation ) is a real constant 


(3) 


which can be taken to be positive, and the negative sign on 
the right is necessitated by the boundary condition at X= so. 


For convenience of discussion one writes (3) as 


Xm — Px (4) 
A } 
(DY) +? aY <0 (5) 
The solution of (4) being obviously 
XK =e>% (6) 


it is that of (5) which is of primary interest. In the follow- 
ing, one will endeavor to study the asymptotic properties of 
any non-trivial solution Y (2, y) of (5) where will 

be assumed to be different from zero, As an apvlication of 
the results obtained in the course of this eG a particular 


problem in atmospheric diffusion will be solved at the end, 
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2, The Asymptotic Behavior of Y (A, y) 


Multiplyine (5) by D , and defining the new variable 4 by 


¥ ie 
NH =f Day (7) 
h 
one has Yt)? sY=0 (8) 
where | 
E(N(¥)) =umpcy) (9) 


In reality, D is different from zero at y = 0 since 
molecular diffusivity is always present, and it must remain 
everywhere finite, Consequently the integral in (7) exists for 
any finite y but increases indefinitely as y-»@. Thus, 
D being a positive anentaey, Y) is a monotonically increasing 
function of y mapping the interval O£y< co into Os <8 ‘ 
Sometimes for convenience a simple functional form is assumed 
for D for instance Dwy™ . But, n_ being usually less 
than 1, the interval (0, co ) is still mapped into itself by 
the transformation from y to r « One will notice in addition 
that since u, ‘OD shane n are all. monotonically increasing 
functions of y ,° g (fy) must be monotonically incressing 
function of 4 . With this in mind, one will consider (8) in 
the interval 0 =1< oO , 

Making the transformations 

Y=GCnH)F(g), &=8&(n) (10) 

and using primes to denote differentiation (with respect to 


y and e ), one has 
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Taking 

eater. (16) 
one has 
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From (13) and (16) one notices that 
» They ekg 
so that € is a monotonic increasing function of r) , becoming 
oO as 4 becomes “, 
It will now be proved that the two terms involving g in 


the parenthesis of (17) vanish asymptotically. One considers 


first the term g’g’g73 . It is sufficient to show that 
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where the lower limit is oo since g"? > 0 as —>), 
Since g is different from zero for all values of 4 different 
from zero, gn2 is finite for such values of " , so that 
the intesral is convergent and s( Fj ) must vanish at cc, 
epee the same argument can be applied to prove that 

g = i Aaa () aS 4 —> 
if m>1l. Now, taking the term g’e-* and differentiating, 
one has 


= eA . 
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The term on the left vanishes asymptotically since 
ES "0 25°05 oO 
and since g is a smooth function. The second term on the 
right has just been shown to vanish asymptotically. Consequently, 
Eu .0 a. A erspoa 
Remembering that O=1F) <o©is mapped into O<E&<~ 
instead of (17) one can study the equation 
FYW+ OX? +9(2)) F =O (19) 
where q(§)-»0 as € ~»0o. It will be assumed that q is 
monotonic asymptotically. After Courant and Hilbert (1931), 
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where & and & are functions of - the asymptotic behaviors of 
which are to be investigated. Calculating F" in two ways 
from (19) and the second equation of (20), one has 
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Rembering that, with q vanishing monotonically for suf- 
ficiently large values of € and with > being positive, the 
quantity q 
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"amplitude" of a solution of (19) is mttun toy at) constant. 
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or if q is of even smaller order, then from (24) it can be 
seen that eae! exists and the distribution of the zeros of F 
will have the even spacing as with respect to oe asymptotically, 
in such a way that after a sufficiently large gero of ° - 3 
all subsequent zeros can be approximatelv located from it by 
using the asymptotic spacing. If q does not satisfy the 
recuirement stated, then the accumulation of the difference 
between < and the actual spacing of zeros will become infinite: 
so that even if the spacing may be approaching ae asymptotically, 
it is impossible to locate from this asymptotic spacing all the 
zeros subsequent to a sufficiently large one, without encounter- 
ing grave errors after sufficiently many locations. What is 
true of the zeros of F is of course also true of those of 
FeA, where A is a fixed number, 

3. %qThe Solution of a Particular Problem in Atmospheric 

Diffusion 
One considers the case where the vertical diffusivity and 

the wind velocity vary as power functions of height, and the 
ground is impervious to vapor, With the vapor concentration 
known to be a certain function of height at x =0 , it is pro- 


posed to calculate the vapor concentration for all positive 


values of x. 
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precisely (u&)2 ° Jes §) f where J, (4 §) are the Bessel 
functions of order +6. Indeed, the principal asymptotic 
properties of the Bessel functions are deduced from those of 
F . The definitions of the Bessel functions are such that %., 
for F Te “e/2" 4 
eft 

The asymptotic properties of F can be utilized, with the 

help of Dirichlet's integral theorem and (33), to furnish in a 


purely formal manner the formula due to McRobert (1931): 
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for Bors and a c(s) vanishing sufficiently rapidly as J7% 
For a rigorous justification of the derivation, however, a few 
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at the moment, Instead, one will proceed with the solution of 
the proposed problem, which will be seen to depend on (35) 
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where and Ag are connected by (31). Let B- we) at 


x= 0. The"density function" (4.) should satisfy 
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In order that the solution be valid, however le(5)| should 
be asymptotically of an order not higher than that of G oo eae ; 
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